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Abstract 

We are concerned with the problem of the decay of a tangle of quantized vortices in He II 
generated by a heat current. Direct application of Vinen's equation yields the temporal scaling of 
vortex line density L ~ t . Schwarz and Rozen [Phys. Rev. Lett. 66, 1898 (1991); Phys. Rev. 
B 44, 7563 (1991)] observed a faster decay followed by a slower decay. More recently, Skrbek and 
collaborators [Phys. Rev. E 67, 047302 (2003)] found an initial transient followed by the same 
classical t -3 / 2 scaling observed in the decay of grid-generated turbulence. We present a simple 
theoretical model which, we argue, contains the essential physical ingredients, and accounts for 
these apparently contradictory results. 
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I. INTRODUCTION AND MOTIVATION 



One of the most important problems of fluid dynamics is the decay of homogeneous 
isotropic turbulence l|, |2j. The problem is particularly interesting in He II (the low temper- 
ature phase of liquid 4 He) because it involves aspects of both classical and quantum physics 
3] . Helium II is a quantum fluid, which can be described phenomenologically, as the mixture 
of two fluid components: the inviscid superfluid of density p s and velocity v s and the viscous 
normal fluid of density p n and velocity v n . Helium's total density is p = p n + p s . Because of 
quantum mechanical constraints, superfluid vorticity is concentrated into thin vortex fila- 
ments of quantized (fixed) circulation n = 9.97 x 10~ 4 cm 2 /s. Superfluid turbulence consists 
of a tangle of such vortex filaments. The intensity of the turbulence is described by the 
vortex line density L (defined as the total vortex length per unit volume); in experiments, 
the vortex line density can be inferred using various methods (e.g. measuring temperature 
gradients, the attenuation of second sound, or the number of ions trapped in the quantized 
vortices). Superfluid turbulence can be easily generated by forcing a sample of helium II 
with moving grids or propellers - techniques which also induce turbulence in classical ordi- 
nary fluids. In this paper we are concerned with thermal counterflow: this is a technique 
without direct analogy in classical fluid mechanics, in which turbulence is caused by the flo! 
w of heat in a two-fluids system. The importance of counterflow turbulence arises from the 
cryogenic engineering applications of liquid helium, because turbulence limits the ability of 
helium II to transport heat. 

The counterflow configuration consists of a channel of cross section A which is closed at 
one end and open to the helium bath at the other end. A heater (electrical resistor) dissipates 
a known power W at the closed end. Normal fluid carries the heat flux Q = W/A = pSTV n 
away from the resistor (where T is the temperature, S the specific entropy and V n the normal 
fluid velocity (averaged over the cross section of the channel); at the same time superfluid 
flows in the other direction to conserve mass, so the total mass flux is zero 

PnV n + PsV s = 0, (1) 

where V s is the average of the superfluid profile. In this way a relative velocity V ns = V n — V s 
between normal fluid and superfluid components is set up which is proportional to the applied 
heat flux, V ns = Q/(p s ST). The pioneering experiments of Vinen [4J determined that, for 
V ns > V\ (where V± is a critical velocity) a tangle of vortex lines of density L = 7 2 fills the 



channel, where the coefficient 7 depends on T. Vinen also wrote a model equation for dL/dt 
which balances growth and decay terms, whose steady state solution (when dL/dt = 0) has 
indeed the observed form L ~ V^ s . The analysis of Schwarz J5], who pioneered numerical 
calculations of superfluid turbulence, confirmed the validity of Vinen's equation, which he 
re-derived from the microscopic laws of vortex dynamics under a number of approximations. 

Further experimental work by Tough jf}] discovered the existence of two regimes of super- 
fluid turbulence: a weak regime (called Tl) for V\ < V ns < V2 characterized by a small value 
of 7, and a strong regime (called T2) for V ns > V2 in which 7 is larger, where V2 is a second 
critical velocity. The nature of the Tl and T2 states is not clear, particularly because the 
transition seems to depend on the channel's geometry; a stability analysis Q| suggests that, 
in cylindrical channels and for V ns > V2, the normal fluid should become turbulent too. 

The decay of counterflow turbulence was experimentally investigated by Schwarz and 
collaborators js, 9], who monitored the vortex line density L as a function of time after the 
heater has been switched off. They discovered that an initial rapid decay of L is followed 
by a slower decay. This result seemed puzzling, as Vinen's equation implies a decay of the 
form L ~ t . In order to make sense of their own experiments, Schwarz and Rozen 8] 
developed a phenomenological model which includes the normal fluid's viscous dissipation, 
and managed to fit the model to the experimental data. 

In the years that followed the attention moved away from thermal counterflow and turned 
to other forms of turbulence. Donnelly and collaborators 10|, 111) found that the decay of the 
vortex line density behind a towed grid follows a L ~ t _3//2 behaviour, which is consistent 
with the classical decay of large turbulent eddies [12] . The result stimulated the study of 
similarities between superfluid turbulence and classical turbulence, and researchers observed 
the same k~ 5 ^ 3 Kolmogorov energy spectra in continuously forced turbulence 13] , the same 
pressure drops along pipes and channels jl^ . and even the same drag crisis for a rapidly 
moving sphere 15) . 

More than ten years after Schwarz and Rozen 8J, Skrbek and collaborators performed 
more experiments on the decay of counterflow turbulence in Prague They found 

that, after an initial transient which seems to depend on the applied heat flux, the vortex 
line density decays as L ~ t" 3 / 2 , the same decay observed in the towed grid experiment, in 



apparent disagreement with Vinen's equation and Schwarz et al. 



9j. The Prague results 



excited the low temperature physics community, because they strongly suggested that the 



3 



decay of counterflow turbulence (which up to that time was generally considered a different, 
non-classical form of turbulence) could be brought into the domain of helium II experiments 



which can be understood in terms o 
^rague classical t~ 3 ^ 2 behaviour 16 



classical fluid dynamics. But the disagreement between 



171 ] and the slower decay observed by Schwarz and Rozen 



8| stood out without any serious attempts to explain it. 

The aim of this paper is to reconcile the experiments of Schwarz and collaborators js, 9| 



with those performed in Prague 



161 ] and show that they are not in contradiction with each 



other. To achieve this aim we shall need to correct the original model of Schwarz and Rozen 
js]. The rest of the article is organized as follows. Section [TT1 introduces and discusses 
Vinen's equation. In Section II III we describe the relevant experiments, while the following 
Section HVl introduces our model, with results given in Section |V] We discuss them and draw 
conclusions in Sectior JVll 



II. VINEN'S EQUATION 



In this section we review the basics of superfluid turbulence in counterflow experiments. 
We make use of the evolution equation for the vortex line density first introduced by Vinen 
4] and then microscopically derived by Schwarz js] using the Local Induction Approximation 
18]. 



In Schwarz's notation, the space curve s = s(£, £) is the position vector along a vortex 
line, where £ is arc length; if we denote the derivative with respect to arc length with a 
prime, then s' is the unit tangent vector, s" is along the normal direction, l/|s"| is the local 
radius of curvature and s' x s" points along the binormal direction. Since the vortices tend 
to grow in the plane perpendicular to the direction of the counterflow velocity v ns = v n — v s 
(Donnelly-Glaberson instability 19j, l20[), the direction which is binormal to the vortices, 
s' x s", tends to be parallel to the direction of the counterflow, v ns . 

According to Schwarz's calculation, the vortex line density equation near equilibrium 
obeys 



dL 
It 



a 



(2) 



where a is a temperature dependent mutual friction coefficient related to Vinen's mutual 



friction coefficient B 



21 



22| by a = Bp n /(2p), 
L- 1 / 2 < s' x s" > and 



< |s"| 2 >= ciL 



2-^1 



(3) 



v n s||v ns is the unit vector and < / >= / f(£)d£/(QL) denotes the line-length weighted 
average of any quantity / where Q is the volume and the integral is along the vortex lines. 
The coefficient (3 is the self-induction parameter, defined by 



K 



/3 = i -ln(L- 1 / 2 /a ), 



(4) 



where a ~ 10~ 8 cm is the vortex core radius. 

A number of authors have dealt with Vinen's equation (j2J) and have attempted to include 
more physical ingredients, taking into account inhomogeneity, anisotropy, rotation and/or 
boundaries. These extensions were made essentially using dimensional analysis, because 
of the difficulty in treating these problems from a microscopical viewpoint 23|, |24j. The 
main physical ingredients which should affect counterflow turbulence are anisotropy and the 
presence of the walls. The former appears in the first term at the RHS of Eq. (jSJ), which, 
according to the Vinen's interpretation, induces growth of the vortex line length. Since this 
paper is concerned with vortex line decay, this term is important only in setting up the 
steady state from which the decay begins, hence it seems sufficient for our purposes without 
any further modifications. In a more exhaustive analysis which is not immediately relevant 
to our aim, the first ! term should be substituted by a tensor, as in Ref 25|, because the 
growth of the vortex line density occurs not only in the direction of the counterflow velocity. 



The channel walls also affect the evolution of the vortex line density 26J]. Mongiovi and 



collaborators 



23 



271 ] extended Vinen's equation (j2J) to include boundaries. The relevant 
dimensionless parameter is the ratio l/D, where i ~ L -1 / 2 is the average inter vortex 



spacing and D is the channel size. They found 23| that for small times the influence of 
the walls is not important and L decays as whereas for long time (in the dilute vortex 
tangle limit) a term which is exponential in time appears. The experimental data which 
we consider {§], are not in this regime: Skrbek 16| at T = 1.6 K, Q/A = 0.08 W/cm 2 
and channel's diameter D = 0.9 cm has i = 7.4 x 10~ 3 cm at t — (see Table HI), and 
£ = 0.031 cm at t = 10 s; Schwarz and Rozen 3 at T = 1.9 K, Q/A = 0.108 W/cm 2 and 
channel's small size D = 1 cm have i = 7.69 x 10 -3 cm at t — (see Table HIT) . and I = 0.3 at 
t = 600 s. We conclude that I < D in the range of interest, hence, for the sake of simplicity, 
we need not to include wall effects in Vinen's equation 

The steady vortex line density solution of Vinen's equation (El) is 



7 V 2 , 



(5) 



T 


Q/A 


V ns (0) 


7 


L(0) xlO 5 


cl 


Ii 


L(20) 


1.6 


0.88 


16.012 


93 


22.176 


0.082 


0.754 


47 


1.6 


0.57 


10.372 


93 


9.3 


0.085 


0.727 


50 


1.6 


0.31 


5.641 


93 


2.75 


0.09 


0.687 


56 


1.6 


0.22 


4.003 


93 


1.386 


0.092 


0.672 


59 


1.6 


0.14 


2.547 


93 


0.561 


0.095 


0.651 


63 


1.6 


0.08 


1.456 


93 


0.183 


0.099 


0.625 


70 



TABLE I: Parameters from the experiments performed by Skrbek et al. LL6J] for different heat 
fluxes: Temperature T (K), heat flux Q/A (W/cm 2 ), steady state counterflow velocity V^ s (0) 
(cm 2 /s), 7 (s/cm 2 ), initial vortex line length £(0) (1/cm 2 ), and parameters cl and In the last 
column we write the value of L at t = 20 sec from our simulations. 



T 


Q/A 


V n8 (0) 


7 


L(0) xlO 5 


CL 


h 


L(1000) 


1.9 


0.065 


0.558 


140.1 


0.061 


0.156 


0.507 


1.17 


1.9 


0.108 


0.928 


140.1 


0.169 


0.15 


0.525 


0.76 


1.9 


0.237 


1.778 


140.1 


0.62 


0.143 


0.551 


1.03 



TABLE II: Parameters from the experiments performed by Schwarz and Rozen 8j at T = 1.9 K. 
The table lists temperature T (K), heat flux Q/A (W/cm 2 ), steady state counterflow velocity 
V ns (0) (cm 2 /s), 7 (s/cm 2 ), initial vortex line length L(0) (1/cm 2 ), and parameters cr and Ii. The 
values of 7 are taken from Ref . 28f] , and the values of cl and Ii are found from Ref . [8| . In the last 
column we write the value of L at t = 1000 sec from our simulations. 

where 7 = Cjj/j3 (with = h/c\) is a coefficient which depends on the temperature, and, 
in some cases, on whether the turbulence is in the so-called Tl or T2 state {(J. We shall 
discuss the issue of 7 in Section IVll 

It is generally assumed that in counterflow turbulence the vortex line density is relatively 
homogeneous, has only a mild anisotropy, is rather featureless (unlike ordinary turbulence, 
where the energy is distributed over the length scales according to the Kolmogorov —5/3 
law), and that the only characteristic length scale of the tangle is l\ this means that both 
the average radius of curvature and the average distance between the vortex lines are of 
the order of magnitude of I. According to Vinen's interpretation, the first term of equation 



6 
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Pn 


Ps 


P 


S 


P 


a 


K 


g/cm 3 


g/cm 3 


g/cm 3 


cm 2 /Ks 2 


g/scm 




1.6 


0.02358 


0.12163 


0.14521 


2824 x 10 3 


1.306 x 10" 5 


0.097 


1.9 


0.06103 


0.08444 


0.14547 


7255 x 10 3 


1.347 x 10" 5 


0.206 



TABLE III: Values of selected parameters in superfluid helium used in our model. 



(J2J) describes to the generation of vortex line due to the relative velocity between normal 
and superfluid component, and the second term models the destruction of vortex lines. 
Essentially, the counterflow supplies energy for the growth of vortex length through the first 
term of Vinen's equation and the tangle converts back its length into energy. In the steady 
state situation, two contributions balance each other, and formula (JSJ) holds true. 

In order to study the decay of the vortex line density in counterflow experiments we need 
dynamical equations for the superfluid helium. The most commonly used model is the two 
fluids model [3|. The two-fluid equations are 



Ps (j^ + v s • Vv s j = -^Vp + p s SVT + F ns , 

fdv n „ \ 

Pn I + V n • VV n I = 

= -^Vp - Ps SVT + pV\ n - F ns , 
P 

V ■ v n = 0, V ■ v s = 0, (6) 

where p is the pressure and p the viscosity. We do not know the precise form of the mutual 
friction force F ns , but, for V ns > Vi, it is reasonable to assume the Gorter-Mellink form: 

F ns ~ p s K«L(v n - V s ) . 



III. EXPERIMENTAL DATA 



Before starting our analysis, we summarize the necessary helium parameters (for p s , p n , 
S and p at the temperature T of a considered experiment see Table IHIj) and the available 
experimental data (Tables HI and HB . 

The Prague experiments were performed in a circular channel (9 cm in diameter and 
13 cm long) at the same temperature T = 1.6 K and at different powers Q/A = 0.88, 0.57, 

7 



0.31, 0.22, 0.14 and 0.08 W/cm 2 . The main finding was that, after a transient t sat , the 
vorticity kL decays as t _3//2 . 

The experiments of Schwarz and Rozen jg] were performed in a rectangular 1.0 x 2.32 cm 
channel which was 24 cm long at T = 1.4 K and Q = 0.043W/cm 2 , T = 1.6 K and 
Q/A = 0.065 W/cm 2 , T = 1.9 K and Q/A = 0.065, 0.108 and 0.237W/cm 2 . Schwarz and 
Rozen presented graphs of the time evolution of the quantity vs t, where L m is L 

times a geometrical coefficient, and /3 is probably taken constant, in the range 0.0 to 100 s, 
except for T = 1.9 K and Q/A = 0.108 mW/cm 2 , for which data are plotted up to 1000 s. In 
this last experiment, Schwarz and Rozen found that the vortex line density initially decays 
as L ~ t -1 , but, after a transient, the decay becomes faster. 

The experiments of Milliken and Schwarz {9] were performed in an open geometry at 
T = 1.33 K, T = 1.45 K and T = 1.60 K. They did not report the heat flux, so we 
cannot model their results. They claimed that the decay of the vortex lines has the form 
L(t) = AkL(0) / (4k + (3L(0)t), where L(0) is the initial value of L, (3 is constant and k = 7.5 
is a coefficient that best fits the experimental data. No data were reported for t > 50 s. 

Not all information which we need in order to set up our model is available in the litera- 
ture: for example, data for L vs t can be read from published figures, but in some cases the 
precise initial value L(0) is not reported, or is not clearly visible in graphs which plot 1/L 



vs t. Moreover, in the Prague experiment 



161 ] the flow channel was cylindrical in shape and 



the planar gold-plated micropore membranes based second sound transducers made the top 
and bottom of another short cylinder oriented perpendicularly to the the channel axis in the 
middle of its length; the channel cross-section was therefore not of exactly uniform shape 
and the measured steady-state vortex line density inferred from the attenuated amplitude 
of standing wave second sound resonance was most likely slightly affected by this geometry. 

To fill such gaps we rely on equation fl5]) for the steady state: given the temperature 
and the initial heat flux, equation (jSJ) enables us to recover the initial value of the vortex 
line density. A problem, however, arises: the parameter^ in equation (JHJ) is not well-known 



(different values are reported in the literature, see [6|, |26|, |28[ and references therein), and 
perhaps is not even unique (in some geometries it depends on the turbulence state being Tl 
or T2). Since our main result does not depend strongly on this choice, we use the values of 7 



recently reported in Ref. 



28 



values are discussed in Ref. 



which we interpret as referring to the Tl state; experimental 
6|, |29( . The choice of 7 is discussed further in SectionlVTl 
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Values of 7 given in Tables |T] and [TT] are therefore taken out from Ref. 28[, and L(0) 
is derived by formula (jSJ). The exact value of (3 can be found from equation (J3J), then 
c l = ftl- The parameters // and Ii/(I\\ — clU) are found using the experimental data for 
— cx,J/) 1//3 respectively from figure 12 and from figure 15 of Schwarz and Rozen's paper 
js]. Parameters J/, In and are defined in Ref. jsj]. The initial counterflow velocity V ns (0) 
is found from the known relation Q/A = p s TSV ns . 

IV. MODEL 

In order to bridge the Prague experiments with those of Schwarz with coworkers, we 
propose the following minimal model, which, we argue, captures the essential physical in- 
gredients. Let x be the direction along the channel and y the direction across it. We start 
from the two-fluid equations ([6]) and the mass conservation equation (CQ). It is fairly certain 
that the normal fluid velocity profile obeys no-slip boundary conditions v n = at the walls 
of the channel, whereas the superfluid velocity can slip along it. 

Our first step is to notice that in a steady state situation, of the order of magnitude, the 
ratio of friction forces and viscous forces acting on the normal fluid is 

friction force napD 2 L , , 

r = 7 « — , 7 

viscous force fi 

where D is the size of the channel. Since r depends only on L, and not on V ns , we can 
estimate r for different temperatures and different vortex line density in the temporal interval 
of interest. At t = 0, the lowest values of L in Tables HI and HT1 are L(0) = 0.185 x 10 5 cm -2 
for Skrbek et al. and L(0) = 0.061 x 10 5 cm -2 for Schwarz et ai, which implies that 
r = 1.6 x 10 4 and r = 3.14 x 10 4 , respectively. At t = 10 s, the Prague experimental data 
show L w 1000 cm -2 , that is r = 873. At t — 600 s in the paper of Schwarz and Rozen, 
L PS (61/3)" 1 , hence r ~ 57. The last columns of Tables HI and HT1 show the values of L at 
larger times resulting from our calculations: we find that the values of r are smaller: r ~ 5 
for Schwarz and r ~ 50 for Skrbek. 

We conclude that, because of the large value of r in the time range of interest (up to 
the time which we shall call t sat ), it is reasonable to assume that the normal fluid profile 



is uniform in y, which is consistent with an experiment of Awschalom et al. 30j. We 



can therefore neglect the relatively thin boundary layer which meets the no-slip boundary 
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conditions at y = and y = D. The same argument applies to the other transverse direction 
z. Thus v n must be independent of y and z. Assuming that the channel is long and thus 
the flow is independent on x, we conclude that v n can be replaced by the scalar quantity V n 
(in the x direction), and that the nonlinear term v • Vv„ at the LHS is zero. This gives us 
an equation for dV n /dt. By applying the same argument to the superfluid equation in ([HD, 
we obtain an equation for dV s /dt. The uniform profiles V n and V s can be identified with the 
cross-channel averaged profiles which satisfy the counterflow condition, thus they are not 
independent: usin! g Eq. ([T|), we reduce the two equations to a single equation for V n , or, 
preferably, for V ns = {p/p s )V n . 

In the experiments, when the heater is switched off, there are short pressure and thermal 
transients during which the driving pressure and temperature gradients change rapidly and 
become negligible; the former is of the order of Atp ~ Yjc\ ~ 10 -4 s where Y is the length of 
the channel and C\ the speed of first sound; the latter is of the order of Atj> ~ Y 2 /x where x 
is the thermal diffusivity, and probably even smaller. We neglect both these transients. We 
also neglect the longer transient during which the vortex tangle depolarizes, and the vortices 
lose the anisotropy imposed by the counterflow: this effect has been already studied 31] and 
involves a correction to L of order unity. Since the anisotropy of the tangle is imposed by 
the counterflow, the corresponding relaxation time is of the same order of that of V ns . 

We conclude that our minimal model reduces to the two equations 

k — aLV ns , (8) 



dt p, 

-<! ( V ns L^ - 

dt V cl 



ah ( V ns L 3/2 - l L A , (9) 



where (3 is given by Eq. (@J, with initial conditions V ns (0) = Q/(p s ST) and L(0) 
(l/( Ps ST)fQ\ 

We stress that our model differs from the model proposed by Schwarz and Rozen 
which consists of three equations for dV n /dt, dV s /dt and dL/dt. Firstly, their model does 
not conserve mass (during the decay of the tangle the counterflow condition still applies, 
because the channel is closed). Secondly, in their model is constant. Thirdly, and more 
importantly, their equation for dV n /dt contains the term —(fi/p n )V n /(D/15) 2 to model the 
effects of viscous dissipation. As we have seen, the viscous forces are negligible compared to 
the mutual friction forces in the temporal range considered in the experiments. Moreover, 
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Schwarz and Rozen openly state that the factor D/15 (rather than, say, D) was chosen 
for the only reason that D/15 gives the best fit to the experimental data. By numerically 
solving their equations, we found that slight changes of this arbitrary factor D/15 produce 
decay curves which are inconsistent with the data. 




FIG. 1: (Color online). Decay of vortex line density L (cm~ 2 ) vs time t (s) for the initial heat flux 
Q/A = 0.88 W/cm 2 at T = 1.6 K, modelling Skrbek's experiment The initial value of the 

vortex line density is 22.176 x 10 5 cm -2 . 



V. RESULTS 

The induction parameter /3 depends on L only via a logarithmic term, so in the first 
approximation it can be considered constant; this is what Schwarz did in his numerical 

n n 

simulations |5| and in the model of Ref. [8[. If /3 is constant, it is apparent that Eqs. (jHJ) and 
([9]) have solutions which scale as L ~ t~ x and V ns ~ r 1/2 . This means that the two terms in 
the bracket of Eq. (Q remain parallel to each other (in log-log axes) and never cross: both 
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FIG. 2: (Color online). Decay of the counterflow velocity V ns (cm/s) vs time t (s) for the initial 
counterflow turbulence intensity corresponding to the same initial conditions of Fig. Q] 

terms scale as t~ 2 . 

However, if L changes by a large amount, which is the case in the experiments under 
consideration, the approximation of constant (3 is not satisfactory. If /3 depends on L as 
in Eq. (j4]), Eq. ([8]) and cannot be solved analytically, so we integrate them numerically 
using the fourth order Runge Kutta method. We find that both L(t) and V ns (t) decrease, 
as shown in Fig. ([1]) and ([2]). 

It is important to notice that, at a certain time t = t sat , the counterflow velocity V ns 
becomes of the order of the typical turbulent superfluid velocity Vt in the vortex tangle, which 
we estimate to be of the order of magnitude of n/i = kL 1 ^ 2 ; hereafter we set Vi = ckL 1 ^ 2 , 
where c is a constant of order unity. Fig. [3] shows the two curves V ns (t) and Vg(t) near 
t = t sat . Clearly for t > t sat the turbulence becomes qualitatively different, because the 
external drive V ns has become as weak as the background turbulence noise. 
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FIG. 3: (Color online). Counterflow velocity V ns (cm/s) (solid green line ) and turbulent velocity 
Ve (cm/s) (dashed red line) vs time t (s), corresponding to the same initial conditions of Fig.[TJ 
For c = 0.7 two curves overlap at t = 0.655 s, which is the value of t sa t plotted in Fig. [4] for 
V ns (0) = 16.0124 cm/s 2 . 



Fig. H] shows our computed values of t sa t for c = 0.6 and c = 0.7 as a function of V^ s (0). 
It is apparent that there is a fairly good agreement with the Prague experiments fit] ]. 

What happens for t > t sat ? When the drive V ns has become very weak, Vinen's equation 
reduces to dL/dt = —((3/cl)L 2 and one would naively infer that L ~ t . A more complex 
scenario is thought to take place, as envisaged by various authors 10L Il2|. which is worth to 
summarise briefly. At the beginning of the decay the superfluid velocity field is essentially 
random, and £ ~ L 1 ! 2 is the only length scale of the tangle. This means that most of the 
kinetic energy is concentrated at large wavenumbers, of the order of As the action of 
the counterflow weakens, it is reasonable to assume that this energy is shifted to smaller 
wavenumbers as in the decay of ordinary turbulence 32] , until the energy- containing eddies 



13 



10 



J 10° ^ 



10' 



: + 








1 1 1 






: o 












- 


" □ 


+ 










- 


















□ 


+ 




O 






- 




□ 


+ 






- 








LJ 


-i- 

5 


n 

I I 












, , i 


+ 





10 



o 



10 1 

V ns (0) 



10' 



?IG. 4: (Color online). Comparison of values of t sa t (s) vs V^ s (0) (cm/s) observed by Skrbek et al. 
lf| (blue crosses) with those predicted by our model with values of t sat for c = 0.7 (red squares ) 
and c = 0.6 (green circles). 

become of the order of the channel's size D. Based on this scenario, the following model for 



the decay of grid turbulence in helium I 
(see also Skrbek, Niemela and Donnelly 



was proposed by Stalp, Skrbek and Donnelly [? 



33| and Skrbek and Stalp 



34|). If U is their speed 

and E their energy, the lifetime of these eddies is r ~ D/U and dE / dt ~ E jr ~ U 3 ~ E 3 ! 2 , 

hence their energy decays as E ~ t~ 2 . If, in analogy with classical turbulence (for which the 

rate of energy dissipation is e = uu 2 where u is the average vorticity and v the kinematic 

viscosity), we identify the average superffuid vorticity with k,L, then e = —dE/dt ~ t~ 3 , so 

u'k 2 L 2 ~ t~ 3 and we conclude that L ~ t~ 3 l 2 . This is the universal decay law, which was 

□ 

first observed by Smith and coworkers [lfj in both towed grid and counterflow experiments 



16 



17J for t > t 



sat- 



and later investigated in detail in Prague experiments 

If we plot the data of Schwarz and Rozen jg] we see in Fig. |5]that, for T = 1.9 K and 
Q/A = 0.108 W/cm 2 , for t > t sat = 600 s the decay is consistent with the same t~ 3 / 2 power 
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* 10' - 
~ 10° t 



FIG. 5: (Color online). Plot of /3L/k vs t (s). The crosses are experimental data from figures 17, 
18, 20 measured by Schwarz and Rozen §| at T = 1.9 K and Q/A = 0.108r W/cmi 2 . The (green) 
line shows the i -3 / 2 dependence. Note the qualitative agreement for t > t sa t = 600 s, as predicted 
by our model. 



law as observed in Prague [16|, |17| , and our calculation predicts t sat = 600 s, provided that 
we take c = 0.77. This agreement between the data reported in Ref. 8] and 16j in terms of 
time dependence has not been noticed before. It is also important to notice that the values 
of t sat are much shorter for the Prague data than for those of Schwarz and Rozen, due to 
the different values of vortex line density in steady state from which the decay originated. 

We have considered the other two cases of Table [Til at T = 1.9 K even though Schwarz 
and Rozen do not plot data for long enough time to reach t sat (up to 100 s). Our model 
predicts t sat = 151.7 s for Q/A = 0.065 W/cm 2 . The important difference (apart from 
the sizes of the channels used, which are, however, nearly the same) between the Prague 
experiments and those of Schwarz and Rozen is thus the vortex line density, which in the 
latter is considerably less, hence the intersection point which defines t sat is not reached. 
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Regarding the experiment at Q/A = 0.237 W/cm 2 , we have not found an intersection 
t sat before t = 1000 s. Probably in this case our model fails because it is too simple: the 
small value L = 1.03 cm -2 in Table HT1 at t = 1000 s implies £ ~ 0.9 cm, which is of the same 
order of the small size of the channel, indeed r w 5, and viscous forces are not negligible in 
this case. 



VI. DISCUSSION AND CONCLUSIONS 

The model which we propose grasps the main behaviour of Skrbek's and Schwarz and 
Rozen's experiments, but the exact value of the ratio c = V ns /nL l l 2 cannot be determined 
for a number of reasons: firstly, our model is too simple to detect fluctuations of the main 
fields; secondly, the completeness of the Hall-Vinen-Bekarevich-Khalatnikov equation is still 



an open question [25|; thirdly, some parameters and initial data are not fully known; fourthly, 
we have neglected any transient in the first part of the decay arising from any diffusion 
process, any decay of normal fluid turbulence, and any possible decay from a T2 state into 
a Tl state. 

Let us consider the state of the superfluid turbulence more carefully. Tough and collabo- 
rators |f| performed many experiments on counterflow in tubes, using different temperatures, 
heat fluxes, tube sizes, tube shapes (circular and rectangular sections) and different mate- 



rials, as described in Ref. 26J and references therein. These studies showed that the values 
of 7 in ([5]) depends not only on T but also on the sizes of the channel and on the heat flux 
Q applied to the sample. In circular tube there are two states (called Tl or T2) charac- 
terised by the value of 7. Since (roughly) 7^2 ~ %Jti, the T2 state corresponds to a larger 
vortex line density (more intense turbulence). The nature of the difference between the Tl 
state and the T2 state is not clear, although it has been suggested [?| that the transition at 
Vns — V2 from Tl to T2 corresponds to a transition to turbulence in the normal fluid (which 
would enhance the intensity of the superfluid vortex tangle). 

Prague experiments were performed in a circular pipe, whose diameter (0.9 cm) was much 
arger than diameters used by Tough and collaborators. If we use Martin and Tough's paper 



261 ] to determinate the critical velocity V2 as a function of D and T, we find V2 = 0.244 cm/s, 
which is smaller than any initial velocity in Prague experiments (see Tabled]). We are thus 
led to suspect that in Prague experiments the turbulence begins the decay from the T2 
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state, although the existing records do not indicate any clear experimental sign that would 
mark the Tl - T2 transition. 



The experiment of Park et al. 35[ used a wide square tube (lxl cm) at relative high 
velocity; they found a single superfluid state, more similar to the Tl state (as confirmed 
by Tough et al. 26[). This result would strongly suggest that in Schwarz and Rozen's 
experiments (rectangular section 1 x 2.32 cm) the decay originated from the Tl state. 

The problem is not only the lack of information about L(0), 7 and the nature of the tangle 
which begins the decay, but also that our model requires values for cl and I;, which, like 
7, are provided only in calculations performed with a uniform normal fluid. In conclusion, 
our model neglects the issue of the Tl or T2 state for lack of precise information (both from 
experiments and from theory) and of consistent numbers to use. 

In conclusion, we have reconsidered the experimental data of the decay of counterflow 
turbulence observed by Schwarz and Rozen [8| and by Skrbek et al. 16j, and compared 
them with our minimal two-fluid model of turbulence decay. The model correctly predicts 
a change in the nature of the vortex tangle at a time t sat which agrees fairly well even 
quantitatively, with the observed onset of semiclassical t _3//2 decay after an initial transient. 
We have found that, contrary to common perception, the early observations of Schwarz and 
Rozen are in qualitative agreement with those of Skrbek et al. Note that the L-dependence 
of f3, ignored by Schwarz and Rozenjs] is essential in our model, because it is responsible for 
the self-induction of the vortex tangle in three dimensions, and hence for the more precise 
description of the decaying tangle. 

Still, our model is perhaps to simple to accurately account for the long time (i.e., t > t sat ) 
asymptotic behavior L oc Dt~ 3 ^ 2 that is displayed over at least an order of magnitude in 



various second sound experiments on decaying counterflow 



10 



] 0. 



lid 



This power law 



behavior seems robust and allowed to determine the values of the effective kinematic viscosity 



^eff 



if the decay law is written in a classical-like form assuming that the Kolmogorov 
K41 form of the energy spectra over length scales up to D (for further details we direct the 
reader to the review of Skrbek and Sreenivasan |40[): 

D(3C) 3 / 2 



L (t) = —± — j=(t + t v 



-3/2 



-t 



-3/2 



(10) 



where C is the Kolmogorov constant and t vo stands for the virtual origin time that in 
most cases can be neglected. Values of is cS extracted this way from Prague experiments 
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on decaying counterflow are in fair agreement with values obtained from the decaying grid 
turbulence in He II. It can be shown that our model would lead to the asymptotic behavior 
L oc Dt~ 3 / 2 if f3 oc y/t/D at late time, or, in terms of vortex line density f3 oc L _1 / 3 .D~ 2 / 3 
rather then j3 oc In (L -1 / 2 /ao) as it was introduced by Schwarz. It is a challenge for the 
future to extend our model (which most likely will have to involve another length scale - an 
outer scale of turbulence D - as in the classical case) in such a way that the experimentally 
observed robust asymptotic behavior will be described more accurately. 
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